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Abstract. This is an investigation of the period- 
luminosity relation of classical Cepheids in samples of dif- 
ferent metallicity. It is based on 481 Cepheids in the Large 
and Small Magellanic Clouds from the blue and red filter 
CCD observations (most similar to Vj & Rj) of the French 
EROS microlensing project. The data-set is complete and 
provides an excellent basis for a differential analysis be- 
tween LMC and SMC. In comparison to previous studies 
of effects on the PL-relation, the EROS data-set offers 
extremely well-sampled light curves and well-filled insta- 
bility strips. This allows reliable separation of Cepheids 
pulsating in the fundamental and the first overtone mode 
and derivation of differential reddening. 

Our main result concerns the determination of dis- 
tances to galaxies which are inferred by using the LMC 
as a base and using two color photometry to establish the 
amount of reddening. We find a zero-point offset between 
SMC and LMC which amounts to a difference between in- 
ferred and true distance modulus of 0.14±0.06 mag in the 
VI C system. The offset is exactly the same in both sets of 
PL-relations — of the fundamental and of the first overtone 
mode Cepheids. No effect is seen on the slopes of the PL- 
relations, although the fundamental and the first overtone 
mode Cepheids have different PL slopes. We attribute the 
color and the zero-point offset to the difference in metal- 
licity between the SMC and LMC Cepheids. A metallicity 
effect of that small magnitude still has important conse- 
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quencies for the inferred Cepheid distances and the deter- 
mination of Hq. When applied to recent estimates based 
on HST Cepheid observations, our metallicity dependence 
makes the low-ifo values (Sandage et al. 1994) higher and 
the high-£zo values (Freedman et al. 1994b) lower, thus 
bringing those discrepant estimates into agreement near 
H ~ 70 km.s" 1 Mpc' 1 . 

Key words: Stars : cepheids - Stars : fundamental pa- 
rameters - galaxies : distances 

* This work is based on observations at the European 
Southern Observatory, La Silla, Chile. 



1. Introduction 

The Cepheid period-luminosity (PL) relation is widely ac- 
cepted as being one of the most accurate primary distance 
indicators to nearby galaxies. It has now been applied as 
far as the Virgo cluster of galaxies (Freedman et al. 1994b, 
Pierce et al. 1994, Sandage et al. 1994, Tanvir et al. 1995) 
and holds the greatest promise to settle the debate over 
the value of the Hubble constant by providing a ±10% ac- 
curate determination (the goal of the HST Key Project). 

The Cepheid variable stars offer a simple way of mea- 
suring distances: their pulsation periods (easy to obtain) 
are strongly correlated with their luminosities; then dis- 
tances are determined from their apparent brightness via 
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the PL relation. The reliability of these distances stems 
from the good theoretical understanding of the Cephcids 
and their PL relation (Iben & Tuggle 1975). Theory pre- 
dicts a small, but not negligible, abundance effect on the 
PL zero point (see Stothers 1988 for review). However the 
ambiguous results of previous empirical tests for this effect 
have led all recent HST studies quoted above to assume 
that the Cepheid PL relation is insensitive to metallicity. 

Three sources contribute to an abundance dependence 
of the PL relation: (1) theory of stellar pulsation, through 
the dependence of period on mass and radius; (2) theory 
of stellar evolution, through the mass-luminosity relation; 
and (3) theory of stellar atmospheres, through line blan- 
keting and backwarming, i.e. the relations between effec- 
tive temperature, absolute magnitudes in bandpasses, and 
bolometric correction. As a result, a metal-poor Cepheid 
is always fainter than a metal-rich Cepheid (at a fixed pe- 
riod and temperature); metal-poor Ccpheids are also hot- 
ter (bluer) on the average. Given the sensitivity of each of 
the above three sources to metallicity, the overall weak- 
ness of the effect on the bolometric PL relation is remark- 
able (Stothers 1988, Chiosi et al. 1993). Theory predicts 
that the slope of the PL relation is nearly independent of 
metallicity, only the zero point is affected. The observed 
PL relations are not bolometric, hence the effect would 
depend on the bandpass. Most of the above predictions 
are valid for a limited range of Cepheid pulsation periods 
— Cepheids with P>50 days are expected (and observed) 
to deviate from a linear PL relation. 

Theory in itself predicts a period-luminosity-color re- 
lation (Sandagc 1958). However we study the PL relation, 
because this is the current one of choice in the determina- 
tion of H . Theory predicts abundance effects on the PL 
relation due to helium (Y) and heavy elements (Z). How- 
ever we only know the heavy elements abundances (Z) in 
Magellanic Clouds Cepheids , hence we speak of metallic- 
ity effects in this paper. 

Previous studies have looked for metallicity effects on 
the Cepheid PL relation since the beginning of the 70s, 
but the issue remains unsettled observationally. Partly to 
blame is the near degeneracy between three important ob- 
served properties of Cepheids: the lines of constant period, 
reddening, and metallicity, in optical and near-infrared 
bandpasses. The lines of constant period represent the 
range of temperatures over which a star can sustain a sta- 
ble Cepheid pulsation in a given mode and period; this 
introduces a natural width (in luminosity) to the PL re- 
lation. The amount of obscuration (reddening) is a more 
serious problem, as it is common practice to derive it from 
the photometry of the Cepheids themselves — thus the 
color difference due to metallicity will affect the reddening 
estimate and the inferred distance (Stothers 1988, Freed- 
man & Madore 1990). 

The first large scale comparison between SMC, LMC, 
and Galaxy Cepheids is due to Paync-Gaposchkin & 
Gaposchkin (1973). Payne-Gaposchkin (1974) concluded 



that there was no evidence for composition differences. 
However, Gascoigne (1974) reinvestigated the apparent 
color differences between LMC and SMC Cepheids as a 
metallicity effect and found that the SMC Cepheids would 
be fainter than LMC Cephcids by 0.1 mag (in V). The 
color shift between LMC and SMC Cepheids was con- 
firmed by Martin, Warren, & Feast (1979) and clearly 
distinguished from differential extinction. Subsequently, 
Iben & Tuggle (1975) and Iben & Renzini (1984) argued 
for a much smaller effect, mostly from theoretical consid- 
erations. Stothers (1988) offered a critical review of all 
these attempts and pointed out the effect of the redden- 
ing correction. In a different approach to the problem, 
Caldwell & Coulson (1985, 1986) merged theoretical and 
empirical relations and individual reddenings from color- 
color diagrams to derive PL relations adjusted for abun- 
dance differences. These fit well their data-set of about 
130 Cepheids in LMC and SMC. Caldwell & Coulson's 
results confirm the existence of a metallicity effect. Their 
very different approach provides no base for a detailed 
comparison with the current extragalactic use (Stothers 
1988, Madore & Freedman 1991). In addition, Caldwell 
& Coulson's sample of about 130 Ccpheids appears to be 
too small for quantifying the tiny effect. This can be seen 
in a recent comparison of primary distance indicators to 
15 galactic and extragalactic objects in which Caldwell 
& Coulson's PL relations {adjusted for metallicity) were 
used (Huterer, Sasselov, & Schechter 1995) — the uncer- 
tainties within and between the four primary indicators 
are larger than the weak metallicity dependence. 

The main reason extragalactic distance measurements 
have assumed that Cepheid luminosity does not depend 
on metallicity is the work by Freedman & Madore (1990). 
They found that any distance differences in three M 31 
fields with varying metallicity are consistent with statisti- 
cal noise, and much smaller than Stothers'(1988) predic- 
tion. The conclusion is based on a sample of 38 Cepheids 
and 152 BVRI observations; the method used is the same 
used in all recent distance determinations. The situation 
was reviewed by Feast (1991), who urged for the need to 
test the metallicity corrections empirically. This was par- 
tially accomplished by Laney & Stobie (1994), who con- 
cluded that their VJHK data on 21 Galactic and 115 
MC Ccpheids do not support the implications of Freed- 
man & Madore's result that intrinsic Cepheid color is in- 
dependent of metallicity. At the same time Gould (1994) 
challenged Freedman & Madore's conclusion by pointing 
out the high degree of correlation among the BVRI mea- 
surements (treated by Freedman & Madore as indepen- 
dent). Gould obtained a PL zero-point shift similar to 
that of Stothers (1988) by reanalyzing the same BVRI 
observations of 36 Ccpheids in M 31. However, he showed 
that the data-set suffers from some systematic uncer- 
tainty, which affects the derived size of the effect. Stiff 
(1995) confirmed independently Gould's conclusions using 
current theoretical evolutionary and atmosphere models. 
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Thus the issue remained unsettled. As a result all recent 
HST measurements of Cepheid distances assume that no 
metallicity effects are present at V and I (Frccdman et 
al. 1994b,Sandage et al. 1994,Tanvir et al. 1995). 

In this paper we use a new data-set of two-color 
photometry (~ 3 million observations) of about 500 
Cepheids in the LMC and SMC to derive the dependence 
of the optical PL relations on metallicity. We find that the 
apparent distance modulus depends on metallicity roughly 
as: 0A([Fe/H}+0.3). The dependence is weaker than some 
previous claims, though by no means negligible. It is de- 
rived from a differential LMC-SMC analysis which incor- 
porates all correlations between Cepheid measurements, 
independently derived metal abundances for Cepheids and 
supergiants, and known limits to the extinction laws. 
The photometric database is a byproduct of the EROS 
microlensing survey (Aubourg et al. 1993a, Beaulieu et 
al. 1995, Beaulieu et al. 1996). 

2. The observations 

The EROS (Experience de Recherche d'Objets Sombres) 
French collaboration (Aubourg et al. 1993a, 1993b, 1995), 
and the MACHO project (Alcock et al. 1993) are both 
searching for baryonic dark matter in the galactic halo 
through microlensing effects on stars of the Magellanic 
Clouds. The EROS CCD equipment has been described 
by Arnaud et al. (1994). The observation, reduction, and 
calibration procedures have been described by Grison et 
al. (1995). All details on the Cepheids discovered and an- 
alyzed in the LMC and SMC are given in Beaulieu et 
al. (1995) and Beaulieu et al. (1996), respectively. 

In brief, all observations were obtained at ESO, La 
Silla using a 0.4-m, f/10 reflector and a 2x8 mosaic of 16 
CCDs. Two broad-band filters were used — Be and Re, 
with central wavelengths which fall roughly between John- 
son B and V and R and /, respectively. The observations 
were obtained during the 1991-92 and 1993-95 seasons. 

Our differential analysis relies on no external zero- 
point, and is thus completed within the Be, Re system. 
However, in applying our metallicity effect to the VI work 
with the HST, we have investigated our transformation 
between Be, Re and V, I. The net transmission of the Be 
band is affected by the blue cutoff of the CCDs sensitivity, 
and is thus much closer to Johnson V than to B. The Re 
bandpass remains intermediate between Cousins R and I. 
The Be filter is broader than Johnson V in a similar way 
as the HST F555W filter is broader than V. The V — I 
and Be — Re colors transform well between each other as a 
result: V-I = 1.02(B E2 -Re2), v = 0.02mags. The LMC 
data from the 1991-92 season was obtain with somewhat 
different set of blue-red filters. Therefore we treated the 
LMC Cepheid data as two different sets and compared 
each of them separately to the SMC data. The 1991-92 
season LMC data is useful because of the large number 
of observations. With thousands of stars available in the 



Cepheids range of magnitudes and colors on the EROS 
CCD frames, constructing an accurate and reliable pho- 
tometric transformation is easy. We transformed all old 
Bei,Rei photometry to the new Bei,Rei system by: 

B E 2-Bei = -0.680 ±0.060- 0.175 ±0.005{B E1 - R E1 ), 

Re2 - Rei = -0.116 ± 0.060 - 0.100 ± 0.005(B S i - R E i) 

The two sets of LMC data are very similar and pro- 
duced the same result in the final differential analysis. 

The identification of Cepheids in both Clouds is based 
on the Fourier components of their light curves and mean 
magnitudes in the expected luminosity range. Several 
Cepheids were excluded from our analysis on the basis 
of strong blending. These could be either actual physi- 
cal binaries, or unresolved stars along the line of sight. 
They are easy to identify by the form of the pulsation 
loop on the color-magnitude diagram (Figure 1) and by 
the accompanying deviation of their mean magnitude and 
color from the rest of the Cepheid sample. Examining the 
pulsation loops shows that background blending for the 
LMC and SMC Cepheids is not a problem (affects less 
than 10% of them), as should be expected for luminous 
stars like them. We investigated several possible sources 



Loops of 1 H Cepheids 




0.4 0.6 0.8 1.0 
B-R 

Fig. 1. The loops which Cepheids complete during a pulsation 
cycle in the temperature-luminosity plane. The form of a loop 
reflects the thermodynamic relation between stellar radius and 
temperature change, plus a small nonlinear effect on the emer- 
gent radiation due to the Cepheid atmosphere. Whenever a 
constant flux is added to that of the Cepheid, the loop looks 
like an "eight" (a bluer blend) or has a "histeresis"-like form (a 
redder blend). One of each examples are shown on this figure 
of first overtone LMC Cepheids — in the upper left and lower 
right sides of the instability strip, respectively. 
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of difference between the LMC and SMC photometry, like 
airmass and the Pinatubo eruption, but found their effect 
to be insignificant. 

3. The concept 

The concept is simple — we have two complete 
Ccpheid samples with known difference in metallicity 
A[Fe/H] LM c-SMC = 0.35 (Spite & Spite 1991; Luck 
& Lambert 1992). This is an average difference between 
spectroscopic abundances of longer-period Cepheids. We 
compare the two samples in the period-magnitude-color 
(PLC) manifold to derive 2 independent sources of differ- 
ence — distance and extinction, and to search for a third 
source — metallicity. 

For each Cepheid we have 3 observed quantities: two 
coordinates and a pulsation period; and 2 observed data- 
sets: complete light curves in Be and Re- The 2 data- 
sets (light curves) give two observed quantities: the time- 
averaged, intensity-weighted mean magnitudes; thus each 
Cepheid has 5 observed quantities in total. The light 
curves also provide independently the means to separate 
Cepheids of different pulsation modes and to treat back- 
ground blending. 

We construct the PL relations in each band, Be, Re, of 
LMC and SMC Cepheids , i.e. the observed mean magni- 
tudes, Qi t k, for each Cepheid of period Pk and each band,i, 
are fit to a linear function: 

Qi,k ~ oti + /3JogP fe , (1) 

where Qj are the zero points and fa are the slopes; within 
our completeness — P < 30days, the PL relations are 
found to be linear. We compare the two Cepheid sets in 
the PL plane. They are offset in luminosity (Qi) due to two 
sources: difference in distance between LMC and SMC, 
and difference in extinction towards LMC and SMC. We 
search for a third source — a term due to the difference 
in metallicity. 

Two PL relations (PL-.Be and PL- .Re) in each of the 
Clouds can give us the distance difference and the redden- 
ing difference, but the metallicity effect remains degen- 
erate. We search to lift this degeneracy by going to the 
pcriod-magnitude-color (PLC) space, and imposing con- 
straints on the interstellar extinction. We fix the mean 
reddening in LMC as E(B-V)=0.10, the value adopted by 
all HST teams, and we adopt independent estimates for 
the foreground reddenings to LMC and SMC. Now, given 
a large observed sample (in a statistical sense), we can 
compare differentially the PLC distributions of LMC and 
SMC Cepheids. Such a comparison will constitute a fit to 
a multi-parameter model. Two of the model parameters 
are metallicity terms. 

We chose as our goal to evaluate the metallicity effects 
on the LMC-based technique used by all HST teams (e.g. 
Madore & Frecdman 1991; hereafter — the modern tech- 
nique) for two reasons. First, the current Ccphcid-based 



H is exclusively derived by it. Second, the EROS Ccpheid 
sample provides an excellent opportunity to accomplish 
that goal, being obtained with two filters (Be, Re) similar 
to the HST F555W, F814W filters, and being LMC-based 
too. 

We emphasize the need to account for the correlations 
between magnitude measurements, which has not been 
done in its applications to-date. This rigorous approach 
was pioneered by Gould (1994); our analysis draws heav- 
ily on it and expands it. Due to the strong correlations 
and near-degeneracies present in the problem, we see no 
alternative to the rigorous approach. 

Hereafter the notation Sx will refer to the metallicity 
dependence of the quantity x in the sense: Axi n f erro d+^ = 
Axtrue, where Ax = x S mc - xlmc- 

4. The method 

Our method applies the same basic techniques for mod- 
eling of data (Press et al. 1994, §15) introduced earlier 
by Gould (1994), but deviates substantially from Gould's 
analysis. The difference is that we solve for a metallicity 
effect which is a function of bandpass/wavelength. There- 
fore we model our data in the three-dimensional PLC man- 
ifold, as opposed to just the PL plane. Thus our analy- 
sis follows the current state of theoretical understanding 
(Stothers 1988; Chiosi et al. 1993; Stiff 1995), and takes 
advantage of the fact that our SMC data-set fills densely 
the PLC manifold. 

Here we present a simultaneous fit of the LMC and 
SMC data. The two (Be, Re) magnitude measurements 
of each Cepheid are treated as correlated with each other, 
but not with the measurements of any other Cepheid . 
Then we can write the x 2 merit function as 

4 N(n) 2 

x 2 = EEE % x h x h (2) 

n— 1 k— 1 — l 

where n = 1,2 correspond to the LMC data (with N(n) 
Cepheids) on the period-magnitude and color-magnitude 
planes, respectively; n — 3,4 correspond similarly to the 
SMC data (with N(n) Cepheids); i = 1,2 are the two 
bandpasses (Be, Re), and where the residuals are defined 
as follows: 

Xl,k = Qi,k - fa + PilogPk), (3) 
x i,k = Qt,k - [at + bi(Q hk ~ Q2,fe)], (4) 
x Ip = Qi,p ~ ( a * + & l °9P P + 7i + l2Ri + ll), (5) 
x i P = Qi,p - [a>i + bi(Qi, P - Q2, P + 72(^2 - Ri) 

+73 1 - 7l) +71+72^+73]: (6) 
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with k,p = l,...,N(n) being the number of Cephcids in 
LMC and SMC, respectively. The covariance matrices of 
the data are then 



N(n) 



n _ \ \ " yn yn 

ij N(n) - 1 ^ ^^j,k 



(7) 



k=i 



for n = 1,2, and 



N(n) 



11 N(n) 
for n = 3, 4, where 



(8) 



x? = 



i 



JV(r. 



(9) 



In equation (2), bij is the inverse covariance matrix of the 
data, b — c _1 . 

The above set of equations illustrates our differential 
analysis of LMC and SMC Cepheids as ensembles of stars. 
Each Cepheid with a period Pfe, and mean magnitudes 
Qi y k is fitted (for each band) to a linear function of the pe- 
riod — the PL relation (pti, Pi), and to the instability strip 
(di,bi). Thus, for the LMC Cepheids we have the resid- 
uals in equations (3) and (4). The residuals of the SMC 
Cepheids will contain four additional linear terms: (1) one 
due to the distance difference (71 = Msmc — Hlmc)] 
(2) one due to extinction — the relative reddening 72 = 
AE(B — V), and Ri being the adopted reddening vector 
(see §6 below); and (3) two (one for each band) due to 
metallicity difference (73). All of these terms simply add 
on to the zero point of the PL relation (for each band), 
and that is what equation (5) is all about. 

The SMC residuals from equation (6) have a slight 
complication, which will become clear here. Let Mj (i = 
1,2) be the absolute magnitudes in Be and Re, re- 
spectively. They define the PL relations for each band: 
Mi(PL) = cti + /3dogP. Then following the metallicity 
term, 73, introduced above, the resulting changes in Mj 
are simply: 8M\ — 73 and <5M 2 = 73. Color is defined as: 
(Qi — ^2)0 = Mi — M2, where the subscript "0" means 
corrected for extinction; hence the color change at a fixed 
period due to the metallicity difference between LMC and 
SMC will be: 



S(Qi -Q 2 )o = (73 1 - 7 3 2 )- 



(10) 



This term, 73 — 73 , will appear in the period-color plane 
together with the reddening (72)-term alone. The equation 
of the linear period-color relation is not added to equations 
(3)- (6) to minimize redundancy; its slope will be simply 

(ft - 

Now we can describe the meaning of equation (6) 
and the residuals in the color-magnitude plane. In the 
color-magnitude plane we have the more complex case, 



where both the color and the magnitude are affected 
by both the reddening and the metallicity. While the 
terms for the magnitude remain the same as in equa- 
tion (5), the terms for the color will come from equa- 
tion (10) for the metallicity, and from the definition 
of extinction color correction. The latter definition is: 
(Qi - 62)0 = (Qi - Ai) - (Q 2 - A 2 ), where the A t are 
the total extinctions in each band, Ai — Eb-vRi, hence 
(Qi - Q 2 )o = (Qi - Q2) - E B _ V (R 1 - R 2 ). Thus equa- 
tion (6) describes the following aspect of our model: more 
extinction makes Cepheids redder and fainter. As far as 
the metallicity terms are concerned, they are not con- 
strained (although theory predicts that more metals make 
Cepheids redder and brighter). We constrain the redden- 
ing term by not allowing unphysical negative values, and 
further by considering the foreground extinction towards 
the SMC as its lower limit (see §5 and 6). 

To obtain the true distance modulus (m — M)o = /i, 
the observed mean magnitude must be first corrected for 
extinction: /i = (<3i-Ai)-Mi = (Q 2 -^2)-M 2 . If A t are 
derived from the observations of the Cephcids themselves, 
the metallicity dependence of the true distance modulus 
will be a combination of the dependencies of the absolute 
magnitudes and the color (Stothers 1988). Using the stan- 
dard definitions for the ratio of total to selective extinc- 
tion, we have: R t = A l /E B -v = (Ri - R 2 )A t / E Qi ^q 2 , 
where E Qi ^q 2 = (Qi - Q 2 ) - (Qi - Q 2 ) . Keeping Ri 
fixed, we obtain for the change of distance modulus with 
metallicity: 



5u 



-6Mi + 



Ri 



Ri — R 2 



-7l + (73 1 -7f; 



Wi-Qj)o 



Ri 



Ri — R2 



(11) 



which applies to the use of a PL relation and is in units 
of stellar magnitude per AZ of metals by mass. Stothers 
(1988) used this case (when Ai are derived from the 
Cepheids themselves) to derive the prototype of equation 
(11) from theoretical considerations. We follow the same 
case (referred to as the modern technique), but derive the 
correction, 5fx, empirically from the differential analysis 
of the magnitudes, colors, and periods of LMC and SMC 
Cepheids. Therefore the application of our method (see 
next section) is restricted to this specific case. We derive 
the metallicity correction as 5fi = Au tIU c — A^i n f orrc d. 

We want to emphasize again that our fit is in the PLC 
space for the sole purpose of finding the residual color 
difference between the LMC and SMC Cepheids ; we are 
not interested in deriving or using a PLC relation. It is 
our way of lifting the degeneracy between the effects of 
reddening and metallicity in the PL planes in our differ- 
ential analysis of LMC and SMC. Freedman & Madore 
(1990) used a different approach to break this degeneracy 
— by setting 71 = 0, i.e. by observing three fields at the 
same distance but different metallicity. They cannot use 
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our method, because their data-set (36 Cepheids ) is too 
small to fill densely the PLC space. In other words, in the 
language of x 2 statistics, their number of data points (36) 
will be comparable to the number of model parameters. 

5. The application 

As described by equations (3)- (6) of the previous section, 
our model of the LMC-SMC comparison has 12 model 
parameters. To solve for all 12 parameters we adopt the 
following external constraints: (1) mean LMC extinction 
E(B-V)=0.10 as used by the HST teams; (2) foreground 
extinctions of 0.06 for LMC and 0.05 for SMC (Bessell 
1991); (3) no depth dispersion in the EROS LMC sam- 
ple; and (4) the line of nodes for SMC from Caldwell & 
Coulson (1986) to derive the EROS SMC sample depth 
dispersion. These constraints are discussed in the next 
two sections. Obviously, there is one more projection — 
on the period-color plane, which is not included in equa- 
tions (3)- (6), but is trivially derived from them. One esti- 
mates the model parameters by minimizing in a statistical 
sense the residuals, given independent constraints on some 
of the parameters. To find the best fit one differentiates 
dx 2 /dA[ = 0, where A[ is the 12 element vector of model 
parameters: 



A = (aii, Pi, ai, b h 71, 72, 73) (i = !> 2 )- 



(12) 



The x 2 merit function of equation (2) can be rewritten in 
a more explicit form for our multidimensional fit: 

4 N(n) 2 12 

* 2 = E E E W"* - E A <-W> (Qi - 

n— 1 k—1 — l 1=1 
12 

x[Ql k - ]T Af^Pk, (Qi - Q 2 ) k )\ (13) 
1=1 

The differentiation yields a matrix equation which can be 
solved for Ar. 



A l = Y J C iq D q , 

q=l 



(14) 



where Ci q is the covariance matrix of the model parame- 
ters vector, 



C-' =B 



4 JV(n) 2 

E/ E/ E/ ^ijfi,k,lfj,k,qi 
n=l k=l 



(15) 



with /™ fe ( being the basis functions from equation (13), 
given explicitly in equations (3)- (6); in a similar fashion, 
the vector D q relates the basis functions to the mean mag- 
nitudes Q" k . 

The procedure requires an iteration (sec Gould 
1994). First, we construct period-magnitude and color- 
magnitude relations for the Cepheids in the LMC and de- 
rive the eight parameters ai,f3i,ai 7 bi from linear fits to 



the data. These are then used to find the covariance ma- 
trices of the LMC data (equation 7). Next, we compute the 
covariance matrices of the SMC data and find the best- 
fit values. We found that the slopes of the PL relations 
do not differ (within the uncertainties) between LMC and 
SMC; we adopt two unique slopes (one for fundamental, 
one for first overtone Cepheids ) for both LMC and SMC. 
At this point we derive A/Xi n f er red, by setting 73 = 73 = 0. 
Then we minimize the residuals along a given reddening 
vector Ry and additionally constraining 72 so that SMC 
Cepheids have no unphysical extinctions (negative). This 
constraint is derived independently of the fit (see next 
section). We use the same form of the Galactic extinction 
law (Cardelli, Clayton, & Mathis 1989) for all samples. We 
iterate this procedure until the best-fit parameters agree 
with the trial parameters used to estimate the covariances. 
Note that the entire procedure is separately done for the 
fundamental mode and first overtone mode Cepheids , and 
thus we have a strong constraint on the final derived pa- 
rameters, which are invariant to the type of Cepheids used 
to derive them. 

For illustration, the covariance matrix of the LMC fun- 
damental Cepheids (equation 7) is: 



0.048 0.041 
0.041 0.037 



(16) 



and for the SMC fundamental Cepheids (equation 8): 



0.103 0.073 
0.073 0.084 



(17) 



The correlation coefficients are obviously very high, as also 
seen in Figs. 5 & 6. 

For Cepheids all at the same distance and selected only 
by magnitude the PL relation is subjected to bias, which 
can be avoided by using the inverse regression (Schechter 
1980, Feast 1995). We do not expect strong magnitude 
bias in our samples, but used both direct and inverse re- 
gressions in our fits. 

The method above contains as a subset the procedures 
which comprise the modern technique for determining 
Cepheid distances (Madore & Freedman 1991). Namely, 
one begins by constructing the PL relations in the LMC, 
then fitting the mean magnitudes of the SMC Cepheids to 
the LMC in each band, and deriving apparent distance 
moduli for each band. These distance moduli are then fit- 
ted to an extinction law with an adopted parameter (Rv) 
to derive the amount of reddening in the SMC data. Fi- 
nally the reddening corrected, true distance modulus is de- 
termined. As brief as it is in this description, the technique 
overlooks several sources of physical scatter and system- 
atics which could be present in the data. We study those 
below, since they are also part of our method to derive 
the metallicity dependence. They provide the independent 
constraints, which make possible the derivation of both 73 
and 73. 
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6. The Extinction Law and Reddening 

In our differential analysis between LMC and SMC the 
effect of interstellar reddening (and obscuration) is repre- 
sented by the term 72 Ri, where 72 = AE(B — V) and the 
reddening vector 

R t = (3.81,2.48). (18) 

We use the independent Galactic extinction law of 
Cardelli, Clayton, & Mathis (1989) and the above values 
correspond to a convolution of the EROS filters transmis- 
sion curves with the Rv = 3.3 extinction law (Fig. 2); here 
R v = A(V)/E(B-V). An extinction law with R v = 3.3 is 

R v = 3.1, 3.3, & 5.0 




0.0 r ................... 1 

1.0 1.5 2.0 2.5 3.0 
Inverse Wavelength (1 /micron) 

Fig. 2. The extinction diagram for three different extinction 
laws — with Rv=3.1 (dashed), Rv=3.3 (solid), and Ry=5.0 
(dotted). Here A(V) corresponds to the amount of extinction 
in the Johnson V band. The two EROS filters are marked on 
the diagram. 

widely accepted for the extinction in LMC (e.g. Freedman 
et al. 1994a), while for our Galaxy Rv = 3.1 is more com- 
mon on the average. There are a few lines-of-sight with i?y 
as high as 5.3, but that's very rare. Physically, high values 
of Ry appear to be related to systematically larger dust 
particles in dense regions. However, an empirical relation 
between the spectral "bump" at 2200 A and Rv seems to 
lack physical grounds (Cardelli, Clayton, & Mathis 1989). 
When used, this relation predicts a high Rv value (as high 
as 5) for the SMC, where the 2200 A feature is known to 
be very weak (Bessell 1991). 

This is a good enough reason for us to investigate how 
our results are affected by the value of the extinction law 
parameter. Therefore we have convolved the transmission 
curves of the EROS filters with R v = 3.1 and R v = 5.0 



extinction laws as well (Fig. 2). After deriving the ap- 
parent distance moduli between SMC and LMC for each 
band and each mode of pulsation, we plot them against 
inverse wavelength at the effective centers of the photo- 
metric bands (Fig. 3). We fit through the points Galac- 
tic extinction laws with different amount of extinction, as 
measured by E(B-V), and with different parameters, Rv- 

E B _ V = 0.08, 0.10, 0.15 



1.4 - 




0.8 - 



0.0 0.5 1.0 1.5 2.0 2.5 
Inverse Wavelength (1 /micron) 



E B _ V =0.10; R v 3.1,3.3,5.0 




0.0 0.5 1.0 1.5 2.0 2.5 
Inverse Wavelength ( 1/ micron) 



Fig. 3. The differential apparent distance moduli for SMC 
Cepheids fitted with Galactic extinction laws: (a) for three dif- 
ferent amounts of extinction, and (b) for three different red- 
dening parameters, Rv- At each band we have two indepen- 
dent values for the distance moduli — for fundamental and 
firts overtone Cepheids , respectively. The error bars on the 
values are as big as the symbols. No correction for metallicity 
effects has been applied at this point. The best fit for SMC 
is E(B-V)=0.09 and 7?y=3.3, but higher values of Rv are not 
excluded. 



Concerning the reddening of the Magellanic Clouds, 
it is important for our study to know the existing ob- 
servational limits for the average values. First of all, we 
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have the limit of foreground reddening. Towards the SMC 
the foreground (due mostly to our Galaxy) is smooth and 
E(B-V) lies between 0.04 and 0.06 mag (Bessell 1991). 
Within the SMC the reddening could vary between 0.06 
and 0.3 mag. On the average the reddening of the LMC is 
very similar: foreground of 0.04 to 0.09, and same within 
(Bessell 1991). These values are derived independent of 
Ccpheid observations and will prove very useful in our 
derivation of a metallicity dependence. Also, when we re- 
fer to Cephcids in LMC and SMC, one should bear in mind 
that the Cepheid samples arc located in the bar of LMC 
and near the central region of SMC. This is particularly 
important in terms of their mean reddenings (as samples); 
the central region of SMC is known to be dusty. 

The extinction diagram in Fig. 3 provides an opportu- 
nity to derive individual reddenings for each Cepheid on 
our list, by fitting an extinction law with a fixed redden- 
ing parameter. The individual Cepheid reddenings are ex- 
cellent for the study of the structure of the PL relation 
(see next section). In deriving individual reddenings we 
required that no Cepheid be assigned a reddening smaller 
than the known foreground reddening to the MC given 
above. As it will become clear from the next section, the 
so derived individual reddenings offer no additional ad- 
vantage in the differential analysis. 

It is well known that interstellar dust absorbs or scat- 
ters light in the optical range, hence extinction is by defi- 
nition a non-negative quantity. In our estimate of the ex- 
tinction, Ai, we need to bear this cutoff in mind, because 
Ai is small towards the LMC and SMC, and for a large 
subsample of our Cepheids we get negative values. Given 
the crucial use of the reddening cutoff in our analysis and 
final fit to the model, it is important to truncate the mea- 
surements of Ai properly, and preserve a useful estimate 
of the measurment uncertainty, a a { in the process of doing 
so. The solution to finding the value of Ai and its error 
together with prior knowledge that Ai cannot be less than 
zero is given by Bayes's theorem. For an estimate of Ai 
with value x and normal error a Xl and our knowledge of 
the distribution of the observed A i7 p(Ai), the probability 
distribution for the true Ai will be given by the Bayesian 
filter: 



p(Ai | x,a x ) = 



pjx | Ai, <7 x )p{Aj) 
p(x) 



p(A t )e 



/ °° p{Ai)e- l ^S^ dA 



A practical, and entirely satisfactory assumption would be 
that p(Ai) is a one-sided gaussian which has a maximum 
at Ai = and declines for large Ai (see Press 1996). 

The conclusion of our analysis is that the reddening 
law in the bar of SMC is not very different from that 
of the LMC with a reddening parameter much closer to 
3.3 than 5. Concerning the absolute value for the LMC 
reddening (currently used: E(B-V)=0.10 and i?y=3.3), 
our mean value, from the individual reddenings, is E(B- 
V)=0.10±0.007. With the derived parameters in the next 



section we could constrain the mean reddening of the SMC 
Cepheids to E(B-V)>0.10. 

7. The PL relations - structure of the scatter 

Now we are ready to apply our method. We start by con- 
structing the PL relations (in each band) for the LMC 
and SMC Cepheids , separately for fundamental and first 
overtone pulsators (Fig. 4). All PL relations are essentially 




Fig. 4. The PL-R E relations for LMC and SMC Cepheids . 

linear in the range of periods in our data. The slopes, (3i 
(eq.l), which are insensitive to moderate amounts of red- 
dening, are the same in LMC and SMC: 
/?i=-2.78±0.16 and /3 2 =-3.02±0.14 for fundamental 
mode Cepheids in LMC, and 

/3i=-2.72±0.07 and /3 2 =-2.96±0.06 for SMC. Similarly, 
for the first overtone Cepheids we have: 
/?i=-3.41±0.22 and /? 2 =-3.43±0.20 for LMC, and 
/?i=-3.46±0.14 and /? 2 =-3.52±0.13 for SMC. This lack 
of metallicity dependence of the slopes, as well as the dif- 
ference between (3\ and /3 2 are in very good agreement 
with theory (Stothers 1988). 

The scatter about the PL relations results from physi- 
cal reasons, as well as statistical noise. On the basis of our 
very well sampled light curves in two bands, we have al- 
ready excluded Cepheids which show evidence of blending, 
as well as stars which are not Cepheids , and can analyse 
the other physical sources of scatter. 

Individually, each Cepheid in the LMC can be off- 
set from the PL relation of the LMC Cepheids due to 3 
independent sources: (1) position in the instability strip; 
(2) differential reddening within the LMC; (3) differen- 
tial distance (depth) within the LMC. These three addi- 
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tional independent sources of scatter contribute to a natu- 
ral width (in L) of the PL relation. They arc manifested in 
the structure of the LMC covariance matrix (equation 7), 
where none of the sources lies in the plane defined by the 
other two sources; hence the covariance matrix would be 
three-dimensional. We illustrate this by plotting against 
each other the magnitude residuals defined in equation 
(3). In Figure 5 we have the diagram of the residuals for 
the LMC fundamental Cephcids . The high level of corre- 
lation is evident. The diagram of the residuals shows that 
reddcing (solid line) within the LMC is not the only source 
of scatter. In Figure 6 we have the same diagram for the 
SMC Cepheids (equation 8); all of the above applies to 
them as well. 



SMC Cep F (R v =3.3; 



LMC Cep F (R v =3.3 



0) 




-1.0 -0.5 0.0 0.5 1.0 1.5 

Res B F 



o 



- K2 



Fig. 5. The residuals in Be and Re, as defined in equation 
(3), plotted against each other for the fundamental Cepheids in 
LMC and their PL relations. The solid line is the reddening line 
for Rv—3.3. The line of constant periods is steeper. Dispersion 
due to distance is even steeper — along the diagonal (dotted 
line) . 



The diagrams of the residuals provide a nice illustra- 
tion to the structure of the PL relations. For example, 
if the PL relation had no dispersion, all points would lie 
in the center (0,0) of the diagram (Fig. 5). Alternatively, 
if there were groupings of Cepheids separated in distance 
along the line of sight, they will appear clearly as clumps 
on the diagram. 

If we deredden each Cepheid as described in the pre- 
vious section §6, we are mostly left with the dispersion 
due to depth (Fig. 7). The reason is that the individual 
dereddening corrects also for the width of the instability 
strip along the lines of constant period, at least partially. 
The correction is partial for two reasons: (1) the redden- 




-1.0 -0.5 0.0 O.i 
Res B E2 



1.0 1. 



Fig. 6. The same as Fig. 5, but for the fundamental Cepheids in 
SMC. The solid line is the reddening line for _fiV=3.3. The 
SMC Cepheids show a larger scatter, which is also closer to 
the diagonal line (due to distance dispersion). 



ing and the P = const, lines are not completely degenerate 
(Fig. 8); and (2) the individual reddening is limited (from 
below) by the amount of foreground reddening (indepen- 
dently known). 

The near-degeneracy of reddening and P = const, lines 
in the color-magnitude plane is worth investigating. We 
find that the slopes of the P = const, lines in LMC and 
SMC are practically the same: 3.3 ± 0.4; while the slopes 
of reddening (B E ) are 2.86 (for i?y=3.3) and 3.53 (for 
i?y=5.0). A limited illustration to this is given in Fig. 8; 
the Cepheids shown have been selected to have exactly the 
same periods (to the second digit) — in practice we use 
the whole data-set to derive the slopes. 

Now we turn to accounting for the depth of LMC and 
SMC. In this we employ for the first time the last two 
observed Cepheid quantities: the coordinates. We plot the 
magnitude residuals versus angle on the sky across the 
SMC (l dc s=0.0389 mag) using the line of nodes for SMC 
from Caldwell & Coulson (1986). They derived them over 
areas which are about 25 times larger than the EROS 
CCD-array areas. The EROS areas cover only the LMC 
bar (which is seen face on) and the SMC bar (seen edge- 

on) and have dimensions of 25 x71 on the sky. That is one 
reason we did not attempt our own derivation of the line 
of nodes of the Clouds. In deriving the depth dispersion 
in the SMC we used Cepheids along the line of nodes and 
orthogonal to it. 

With the width of the instability strip derived from the 
color-magnitude plane to be ±0.22 mags, we derive the 
depth dispersion in our SMC sample to be ±0.10 mags. 
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Fig. 7. The same as Fig. 5, except that each Cepheid has been 
dereddened individually using the Galactic extinction law with 
_Ry=3.3. In comparison to Fig. 5, the dispersion has decreased 
and has shifted slope to the diagonal. All these Cepheids are 
in the bar of the LMC, which is face on, and that explains the 
very small depth dispersion we see. Conversely, in SMC we see 
a substantial depth (Fig. 6). 



P = const. 




0.50 -0.36 -0.22 -0.08 0.06 0.20 

Fig. 8. An excerpt of the Cepheid instability strip showing 
groups of Cepheids which happen to have the same periods — 
from top to bottom: 7.5, 5.5, 4.2, and 3.1 days, respectively. 
The open symbols correspond to LMC Cepheids , the rest — 
to SMC ones. The Cepheids shown have been dereddened as 
an ensemble in each MC, with Rv~S.3. The reddening lines 
for two values of Rv are shown in the lower left corner. 



Our uncertainty in the value of Rv limits our ability to de- 
rive the internal amount of reddening within SMC to the 
ranges already known (§6). However, it provides a very 
valuable constraint from below on the overall reddening 
towards SMC as E(B-V)>0.10±0.01; otherwise half of the 
Cepheids will have negative extinctions, which is unphys- 
ical. We use this value to constrain j 2 m our multidimen- 
sional minimization (equations 12-15). 

Finally, we would like to point out that the EROS 
sample of SMC Cepheids happen to be mostly from the 
far arm and the main bar (see Figs.7&9 of Caldwell & 
Coulson 1986). Hence the sample is not likely to define the 
centroid of SMC well, and would result in a larger distance 
modulus difference with the LMC (by about 0.15 mags). 
This is of no consequence to our differential analysis here. 

The diagrams of the residuals are independent of dif- 
ferential reddening or metallicity between the MCs, how- 
ever understanding the structure of each PL relation and 
quantifying the sources of scatter improve the reliability 
of the differential analysis. When we compare the LMC 
and SMC PL relations, the above three sources of scatter 
are not involved (as long as we have a statistically large 
sample of data points, i.e. Cepheids) — we use ensemble 
quantities. Then it is only consistent to use the "ensem- 
ble" reddening for each MC, the individual reddenings will 
introduce a bias. Other sources of scatter have been pro- 
posed in a critique of the PLC relation — mass loss and dif- 
ferent strip crossings (Stift 1982, 1995). We do not identify 
another source in our data and in the remaining scatter; 
such effects must be subtle for the PL relation alone. 



8. Results 

8.1. The effect of metallicity 

We apply the procedure in § 5 to the data-sets of funda- 
mental mode and first overtone mode Cepheids separately 
and obtain very similar results. There are enough first 
overtone Cepheids , 27 and 141 respectively, to make this 
result statistically significant. Thus, the two sets of PL re- 
lations give the same distance modulus difference (SMC- 
LMC) of 0.86±0.045 mag (before corrections for reddening 
and metallicity) . This illustrates for the first time the use- 
fulness of the PL relation of first overtone Cepheids for 
determining distances. 

The model parameters derived from the fits of the 
different datasets are as follows: a^,/3j= 17.61±0.035,- 
2.72±0.07; 17.74±0.029,-2.95±0.06; a' t ,h= 14.49±0.13, - 
11.78±0.94; 14.36±0.12, -12.77±0.90; 71= 0.62±0.04, 72 = 
0.010±0.01, 0.06±0.01, and 7 |= -0.01±0.01. This is 
the fit with the LMC 1993-94 data; the zero points a-, 
a\ refer to the SMC linear fits. The correction due to 
the metallicity dependence of the inferred distance mod- 
ulus by the described technique is then Sfj, — A/j, true — 
Afx mferred = -0.145 ± 0.06mag. If the LMC 1991-92 
dataset is used, the color shift is smaller by 0.02 mags 
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to 73 —73= 0.05±0.01, but a corresponding change in the 
PL relations leads to 71= 0.63±0.04, 72= 0.025±0.015, 
7 3 1 =0.05±0.01, and 7 f= -0.01±0.01. Therefore, 8fjt = 
—0.139 ± 0.06mag. In the error budget, in order of im- 
portance, the sources are: (1) the photometric transfor- 
mation between the 1991-92 and present data (0.06 mag); 
(2) the scatter in the PL relations (0.04 mag); (3) the un- 
certainty in the reddening (and depth) estimate in SMC 
(0.03 mag). The first error source is not applicable to the 
LMC 1993-94 dataset. 

The above result is for a fixed reddening vector, Ri 
(equation 18). If we change the parameter to i?y=5.0, i.e. 
i£,=(5.51, 3.95), the above result for the metallicity effect 
is no longer unambiguous. The reason for this is that now 
it is possible to satisfy the constraint from the foreground 
SMC reddening and correcting for the offset of the SMC 
Ccpheids in the color-magnitude plane. In other words, if 
in SMC R v =5.0 and E(B-V)=0.04, one need not invoke 
a metallicity effect in the differential comparison to LMC. 
The fact that the SMC Ccpheids appear bluer than the 
LMC Cepheids is then explained with the small amount 
of extinction and steep extinction law. 

While a value of Ry >5.0 is not excluded (but un- 
likely) by our data, an extinction of E(B-V)=0.04 is very 
difficult to accommodate in our analysis of the magnitude 
residuals (see §7). Our mean value is 0.125±0.009, and 
that is in very good agreement with a number of indepen- 
dent estimates, as reviewed by Bessell (1991). 

Therefore we consider this alternative explanation 
(very high Ry and very low reddening) as very unlikely. 
We propose instead that the observed effect is due to the 
known difference in metallicity between the Cepheid sam- 
ples of the LMC and SMC: 0.14 mags for a factor of 2 
difference in heavy elements. The sign of the effect would 
lead to the distance of a metal-poor sample to be overes- 
timated. 

In applying the above result to other galaxies, we need 
to address two issues: (1) about the form of the metal- 
licity dependence of inferred Cepheid distances, and (2) 
about the role of the helium abundance (Y). The form 
of the dependence can be inferred from current theoret- 
ical models. Photometric colors and bolometric correc- 
tions (from model atmospheres) have linear dependen- 
cies on log( q 16 ), where Z is the abundance of metals 
(by mass) and Z=0.016 in the Galactic Cepheids. Mass- 
luminosity relations, hence — Mi, seem to depend lin- 
early on Z (Chiosi, Wood, & Capitanio 1993; Buchler 
et al. 1996). In addition, these models show a significant 
linear dependence on the abundance of helium (Y). We 
have currently no data on AY between the LMC and 
SMC Cepheids, therefore in applying our metallicity de- 
pendence on distance to other galaxies, we have to as- 
sume that (AY/AZ)smc-lmc = (AY/ AZ) galaxy _ LM c- 
That assumption is strongly supported by the current level 
of understanding and measurement in stellar populations 
relevant to our study (e.g. Olive, Skillman, & Steigman 
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1996), and will not affect any conclusions with regards to 
distances. We also want to point out that the above needs 
to be born in mind in comparisons to results from theo- 
retical modelling, e.g. Chiosi et al. (1993) use Y=0.27 for 
both LMC (Z=0.008) and SMC (Z=0.004) in Table 15, 
which may not be a good assumption. 

In view of the above and our model, we de- 
rive the following metallicity dependence of inferred 
Cepheid distances: 

5p = (0.441°;*) log-^—, (19) 

^LMC 

where Z is the abundance of metals (by mass) in the stud- 
ied Cepheids and Zlmc=0-0085. The metallicity depen- 
dence is valid in the spectral region covered by the EROS 
filters. Due to these filters transformation properties into 
the standard BVRI system, and the application at hand, 
we have derived it for V, and V-I use in particular, where 
V-I = \m(B E - Re),& = OMmags. The above metal- 
licity dependence applies to distances inferred by differ- 
encing against the LMC Cepheids and Cepheid derived 
reddenings (what we refer to as the modern technique), 
and R V (LMC)=3.3. 

The metallicity dependence we find is two times less 
steep than that found by Gould (1994) from Freedman & 
Madore's M31 data (0.88±0.16 [Fe/H]). It is also smaller 
than Stothers' prediction (28.7 SZ, for 5Y/SZ = 3.5). 
However our result is in perfect agreement with the metal- 
licity effect seen in Galactic Cepheids as a function of 
galactocentric distance. The observed abundance gradi- 
ent for Cepheids is S[Fe/H]/SR GC = -0.07 ± 0.02fcp C - 1 
(Giridhar 1986). Two independent Cepheid samples show 
the same effect of progressively bluer color with increas- 
ing galactocentric distance (Caldwell & Coulson 1987, and 
Gicrcn et al. 1993), which converts to the following values 
for 73- - 7I (of equation 10) in units of [Fe/H]: 0.29 ± 0.05 
for (B-V), and 0.20 ± 0.05 for (V-I). For (V-I) of our sys- 
tem we obtained 0.20 ± 0.02. Exactly the same value is 
obtained (and a nice illustration of the effect) from the 
compilation of LMC and SMC Ccpheids in the study by 
Di Benedetto (1995) - see his Fig.l & 2, and also 8. 

Our comparison to the work of Caldwell & Coulson 
(1986) in BVI and of Laney & Stobie (1994) in VJHK 
is limited for at least two reasons. First, we cannot repro- 
duce their reddening determinations with only two bands. 
Second, we are unable to comment on the theoretical mod- 
els which they used as constraints. Nevertheless, we find 
a similar intrinsic difference in V — I colors between LMC 
and SMC as Caldwell & Coulson; our value is two times 
larger. We also agree with Laney & Stobie's main conclu- 
sions, despite the uncertainty of a J to 7 transformation 
needed for such a comparison. 

The metallicity dependence is not strong enough to 
upset the overall agreement of the four primary distance 
indicators in fifteen local distances (Huterer, Sasselov, & 
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Schechter 1995). Similarly, the discrepant distance mod- 
uli involving IC 1613 noted by Gould (1994) are less dis- 
crepant with our weaker metallicity dependence; the RR 
Lyrae distances may be partly responsible for the remain- 
ing difference. 

8. 2. Implications for H 

Despite its relative weakness, this metallicity dependence 
has a significant effect on extragalactic distance measure- 
ments. The effect of metallicity on Cepheid luminosity, 
when not accounted for, can lead to two types of errors in 
distance measurements (Gould 1994): incorrect Ho value, 
or a dependence of Ho on distance appearing as proper 
motions or Malmquist bias. Here we illustrate the effect 
of our metallicity dependence on the value of H , as based 
on Cepheids . 

Recent efforts to determine H from Cepheid distances 
have focused on: (1) the Virgo and Fornax cluster galaxies 
(Frcedman et al. 1994b, 1996), (2) the Leo I group, con- 
taining ellipticals (Tanvir et al. 1995); and (3) the parent 
galaxies of supernovae of type la (Sandage et al. 1994). 
All these studies use the same modern technique with the 
LMC as a base and all the same initial assumptions; we 
share them in our analysis. Yet they result in three differ- 
ent values of Ho, with hardly overlaping error bars. The 
abundances of metals in the young populations of the host 
galaxies for these studies differ and, apparently, in a sys- 
tematic way. Therefore we propose that the metallicity 
dependence may be responsible for most of this discrep- 
ancy. To illustrate our idea, below we apply Equation (19); 
please note that the uncertainties in the mctallicities of 
the extragalactic Cepheids are not yet well known, hence 
we do not try to provide rigorous error estimates to the 
corrected values. 

With the first approach and HST VI photometry 
of Cepheids in Ml 00, Frcedman et al. (1994b) derive 
Ho = 80 ± 17 km.s" 1 Mpc^ 1 . For the metallicity of the 
Cepheids we adopt [Fe/H]=+0.1 (Z=0.021), derived from 
the abundances of H II regions (Zaritsky et al. 1994), as- 
suming a solar [O/Fe] ratio. With the uncertainty range 
for the slope of our metallicity effect, this leads to H a = 
76 — 70 km.s -1 Mpc^ 1 . With some of the color difference 
due to metallicity effects, the estimate of interstellar red- 
dening will also change (the metal-rich M100 Cepheids 
are intrinsically "redder" than the LMC ones) — e.g. from 
E(V-I)=0.13±0.06 (after Ferrarese et al. 1996) we obtain 
0.05 mags. Recent interim results reported on half-dozen 
more galaxies (Freedman, Madore, & Kennicutt 1996) in- 
clude mostly spirals which are not as metal rich as M100, 
e.g. M101 has LMC abundances in its outer field. All this 
brings their interim result down to 73±10 from their initial 
MIOO-based value, as we would expect from our metallic- 
ity effect. A complete analysis should be done when all 
data becomes available. 



With the second approach and HST VI photometry 
of Cepheids in M96, Tanvir et al. (1995) derive Hq = 
69 ± 8 km.s -1 Mpc -1 . The abundance estimate for M96 
comes form Oey & Kennicutt (1993), slightly below so- 
lar at [Fe/H]=— 0.02. This implies a small change in the 
Cepheid distance to M96 and a decrease in the value of 
H by 3—5 km.s -1 Mpc^ 1 . The reddening changes from 
E(V-I)=0.09±0.10 to 0.03 mags. 

With the third approach and HST VI photome- 
try of Cepheids in IC4182 and NGC5253, Sandage et 
al. (1994) derive H = 55 ± 8 km.s" 1 Mpc- 1 . Abun- 
dances in NGC5253 have been measured and discussed 
by Pagel et al. (1992); discussion of abundances in IC4182 
is given by Saha et al. (1994) - wc adopt [Fe/H]=-1.3 
and caution on the large uncertainties. With this metal- 
licity the above estimate of H should be increased 
by 6 to 13 km.s -1 Mpc- 1 . The reddening estimates 
change from E(V-I) = -0.11±0.07 and 0.03±0.15, to E(V- 
I)=0.09±0.07 and 0.23±0.15, respectively (these metal- 
poor Cepheids are "bluer" than the LMC ones). Here we 
should point out that the interpretation of the light curves 
of SN la themselves has been improved considerably by 
Ricss, Press, & Kirshner (1995), leading to an additional 
systematic shift upwards by 6 km.s -1 Mpc- 1 . This brings 
about the estimate to H ~ 70, which is also in good agree- 
ment with the theoretical calibration of SN la and the 
value H = 67 ± 9 km.s" 1 Mpc' 1 by Hofiich & Khokhlov 
(1996). The new value of H = 58 ± 4 km.s" 1 Mpc- 1 
by Sandage et al. (1996) includes two SNe la in spiral 
galaxies, in which Cepheid metallicity is most likely sim- 
ilar to that of the LMC. Note, however, that the two 
SNe - 1981B (NGC4536) and 1990N (NGC4639) are by 
0.2—0.4 mags dimmer than the three SNs in IC4182 and 
NGC5253. Therefore, the result we obtained above by ac- 
counting the effect of metallicity on the Cepheid distances 
for each galaxy, remains virtually unchanged, and so does 
the consistency with the results of Riess, Press, & Kirsh- 
ner (1995). 

The metallicity dependence we found from the 
LMC/SMC analysis brings all the derivations of H to 
good agreement. It also leads to reasonable amounts of 
reddening, especially regarding the negative values ob- 
tained otherwise in metal-poor samples. 

Secondary distance indicators, like the Tully-Fishcr 
(TF) relation, the Surface Brightness Fluctuations (SBF) 
method, and the Planetary Nebulae Luminosity Func- 
tion (PNLF) method, have their zero points calibrated 
with Cepheids. The TF relation is calibrated primarily 
in M81 (Freedman et al. 1994a) with 25 Cepheids. The 
metallicity of these Cepheids is high, [Fe/H]«0.05 (Gar- 
nctt & Shields 1987; Zaritsky et al. 1994). The SBF and 
PNLF zero points come from M31, and in a similar fashion 
— from metal-rich Cepheids (Freedman & Madore 1990). 
Given our metallicity effect, we are not surprised to note 
that all these methods give a high value for H - 
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9. Conclusions 

We use the new EROS microlensing survey data-set of 3 
million two-color observations of about 500 Cepheids in 
the LMC and SMC to search for and derive the depen- 
dence of the optical PL relations on metallicity. We find 
that: 

(1) The PL relations for both types of Cepheids have the 
same zero-point offset (but no slope difference), which is 
attributed to the effect of metallicity under a reasonable 
assumption about the extinction law. As expected from 
theory, this effect of metal content is manifested in a color 
shift of the instability strip. It amounts to about \th of the 
strip width, and we could detect it unambiguously thanks 
to the large number of Cepheids and extremely well sam- 
pled light curves. 

(2) With the known ensemble difference in metal content 
between LMC and SMC Cepheids , we derive a linear rela- 
tion between the distance modulus correction and metal- 
licity: 

5/i = (0.44lg;J) log-?— , 

^LMC 

It applies to distances which are inferred by using LMC as 
a base and using two color VI photometry of the Cepheids 
to establish the reddening. The linearity of metallicity de- 
pendence is a good assumption, but needs to be confirmed 
empirically outside the range of application (a factor of 
few lower than SMC and higher than the Galaxy). 

(3) The first overtone Cepheids have PL relations which 
provide distances fully consistent with the PL relations of 
fundamental mode Cepheids . 

We use two color bands closely spaced in wavelength 
because of availability, not their desirability; ideally one 
would like to use at least three bands, one of them in the 
near-infrared. 

Our result can be applied to the long-standing discrep- 
ancy between the low-i/o scale and the high-i?o scale. The 
host galaxies on which each of these scales relies appear 
to have systematically different metallicities. A simple ap- 
plication of our correction to several recent derivations 
makes the \ow-H values (Sandage et al. 1994) higher and 
the high-i?o values (Freedman et al. 1994b) lower, thus 
bringing those discrepant estimates into agreement near 
H a ~ 70 km.s^ 1 Mpc- 1 . 
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